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Position and energy-resolved particle detection using phonon-mediated 
microwave kinetic inductance detectors 

D.C. Moore,^'^) S.R. Golwala/ B. Bumble,^ B. Cornell/ P.K. Day,^ H.G. LeDuc,^ and J. Zmuidzinas^'^ 

^■'Division of Physics, Mathematics & Astronomy, California Institute of Technology, Pasadena, 

CA 91125 

^^ Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109 

We demonstrate position and energy-resolved phonon-mediated detection of particle interactions in a silicon 
substrate instrumented with an array of microwave kinetic inductance detectors (MKIDs). The relative 
magnitude and delay of the signal received in each sensor allows the location of the interaction to be determined 
with <1 mm resolution at 30 keV. Using this position information, variations in the detector response with 
position can be removed, and an energy resolution of cfe = 0.55 keV at 30 keV was measured. Since MKIDs 
can be fabricated from a single deposited film and are naturally multiplexed in the frequency domain, this 
technology can be extended to provide highly-pixelized athermal phonon sensors for ~1 kg scale detector 
elements. Such high-resolution, massive particle detectors would be applicable to rare-event searches such as 
the direct detection of dark matter, neutrinoless double-beta decay, or coherent neutrino-nucleus scattering. 



Next generation rare-event searches such as the di- 
rect detection of dark matter require large target masses 
(~10^ kg) with sub-keV energy resolution. This re- 
quires increasing the mass of current solid-state, cryo- 
genic experiments^'^ by 2 orders of magnitude, while 
maintaining the background-free operation of existing de- 
tectors. Reducing the cost and time needed to fabricate 
and test each detector element is necessary for such large 
cryogenic experiments to be feasible. 

Detectors that measure both the athermal phonons 
and ionization created by a particle interaction have 
demonstrated sufficient background rejection to enable 
next-generation experiments^. Microwave kinetic induc- 
tance detectors (MKIDs)^'^ offer several advantages for 
providing athermal phonon sensors in large experiments 
relative to the transition edge sensor (TES)-based de- 
signs currently in use^'^'^. MKIDs can be patterned from 
a single deposited aluminum film, with large (>10 jim) 
features, significantly reducing fabrication time and com- 
plexity. Since MKIDs are naturally multiplexed in the 
frequency domain, hundreds of sensors can be read out on 
a single coaxial cable, enabling a more granular phonon 
sensor that is expected to provide enhanced background 
rejection. In addition to dark matter direct detection, 
high-resolution, massive particle detectors are applicable 
to the detection of neutrinoless double-beta decay^ and 
coherent neutrino-nucleus scattering^. 

Previous designs^^^ attempted to absorb the incident 
energy in large-area collectors coupled to smaller vol- 
ume, distributed MKIDs. Although separating the ab- 
sorber and sensor allowed increased sensitivity by con- 
centrating the absorbed energy, test devices suffered from 
poor transmission of quasiparticles from the absorber 
to sensor. Here we present a simplified design that 
eliminates the absorber by directly collecting the en- 
ergy using large- area MKIDs. A similar design devel- 
oped independently by Swenson et al.^^ has been used 
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to demonstrate time-resolved phonon-mediated detection 
of high-energy interactions from cosmic rays and nat- 
ural radioactivity using MKIDs. We have previously 
demonstrated energy-resolved detection of 30 keV x- 
rays^^ using phonon-mediated MKIDs fabricated from 
high-resistivity nitride films^^, but the energy resolution 
was limited to a^; = 1.2 keV by systematic variations 
in the reconstructed energy with interaction location. In 
this work, we present results for devices fabricated from 
aluminum films that exhibit improved uniformity. These 
devices allow the location of the interaction to be de- 
termined and a position-based correction to be applied, 
recovering a reconstructed energy resolution within 40% 
of the baseline resolution. 

To calculate the expected detector resolution, we as- 
sume that the phonons uniformly illuminate the sub- 
strate surface and that incident energy is converted into 
quasiparticles in the inductive section of the resonators 
with overall efficiency rfph. In this case, the phonon life- 
time, Tphi is determined by the substrate geometry and 
fraction, rjfuu of the total surface area, Aguhi that is cov- 
ered by MKIDs. To resolve the rising edge of the phonon 
signal requires that the quality factor Q < 10^ for res- 
onant frequencies /o ~ 3 GHz to give a resonator re- 
sponse time Tres = Q/'t^Iq ^ 1 MS. For rjfm ^ 5% as 
in existing TES-based designs^, Tph ^ 1 nis, while quasi- 
particle lifetimes, r^p, in our films are measured to be 
10-100 /is. Given these characteristic times, we assume 
Tres <^ Tqp <C Tp^. The rcsouator quality factor, Q~^ = 
Q~^ + (5~^, depends on both the coupling quality factor, 
Qc, due to energy losses through the coupling capacitor 
and the internal quality factor, Q^, due to all other losses 
(e.g., dissipation due to quasiparticles). For our devices, 
typically Qi > 10^, so Qc <C Qi is required for sufficient 
sensor bandwidth to resolve the phonon rise time. We as- 
sume the steady-state quasiparticle number is dominated 
by the readout power ^^, Preadi (i-e., thermal quasiparti- 
cles or those generated by stray pair-breaking radiation 
are negligible) so Nqp = rjreadPreadTqp/A, where A is 
the superconducting gap and rjread is the efficiency for 
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FIG. 1. (a) Array layout. The resonators are arranged in an 
offset grid and coupled to a coplanar strip feedline. The num- 
bers denote the ordering of resonators by frequency, which 
were evenly spaced by 10 MHz so that the entire array fit 
within 200 MHz of bandwidth around 4.8 GHz. The interac- 
tion locations for the event from Fig. 2 determined from the 
energy partitioning (x) and timing delays (O) are also shown, 
(b) Geometry for a single resonator. 



the readout power to create quasiparticles in the res- 
onator. While the detailed mechanism for quasiparticle 
generation from the sub-gap readout phonons remains 
unclear^^, in the most conservative case we can assume 
Vread = 1- Finally, we assume that frequency-based read- 
out is used with the noise dominated by the HEMT am- 
plifier with noise temperature, Tn- 

Given these assumptions, the expected energy resolu- 
tion for an MKID-based phonon- mediated detector is^'^^: 
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where pt is the probability for a phonon to be transmit- 
ted from the substrate to the MKID, a is the fraction of 
the total inductance of the resonators due to kinetic in- 
ductance, Ap5 is the mean free path for phonons to break 
Cooper pairs in the MKID^^'^^, and A^o is the single-spin 
density of states. The ratio of the frequency to dissipa- 
tion response is defined as /3, while Si is a dimensionless 
factor of order unity describing the dissipation response 
at temperature, T^^'^^. 

Figure 1 shows the design of prototype arrays consist- 
ing of 20 lumped-element resonators^^, each with a 70 /am. 
wide meandered inductor (with total area Aind = 1 mm^) 
connected to an interdigitated capacitor with 20 /im spac- 
ing between fingers. Devices were fabricated on 1 mm 
thick, high-resistivity Si substrates (p > 5 kQ cm), which 
were deglazed in dilute hydrofluoric acid (HF) just prior 
to vacuum processing. The resonators were patterned 
by inductively coupled plasma etching (BCI3/CI2) of a 
single, 25 nm thick aluminum film deposited by dc mag- 
netron sputtering at ambient temperature. The sub- 
strates were mounted to a testing enclosure with GE var- 
nish at the corners of the chip, with total contact area 
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FIG. 2. (Color online) Pulses observed in each resonator for 
a 200 keV energy deposition. The pulse heights, giving the 
phase shift of the signal transmitted past the resonator in de- 
grees, have been calibrated to remove responsivity variations, 
and a 200 kHz low-pass filter has been applied. The numbers 
denote the corresponding positions shown in Fig. 1 for the 
resonators closest to the interaction site, which have large, 
prompt responses. The resonators farther from the interac- 
tion have smaller, delayed responses, (inset) Optimal filter fit 
(solid) to the pulse in the primary channel. The pulse contains 
a prompt component (dash-dotted) that decays away with the 
quasiparticle lifetime as well as a delayed component (dashed) 
that decays away with the phonon lifetime. 



of 8 mm^. An ^^^I source was mounted to illuminate the 
substrate face opposite the resonators, and the enclosure 
was cooled to 50 mK in a dilution refrigerator. 

The measured Qc varied from 5x10^-5x10^, leading 
to a corresponding variation in the responsivity. The 
relative responsivity was calibrated using ^200 events 
from cosmic rays or natural radioactivity that deposited 
>500 keV in the substrate. An example 200 keV inter- 
action is shown in Fig. 2. At times longer than ~100 /is 
after the interaction, the initial position-dependent infor- 
mation has decayed away and the phonons have diffused 
sufficiently that they are uniformly distributed through- 
out the substrate. The phonon velocity is 6-9 mm//is, 
depending on polarization^^, indicating that the average 
phonon has traversed the substrate more than 20 times in 
this period, and all position-dependent information has 
been lost. By matching the amplitude of the pulses' ex- 
ponential tails at long times, the relative responsivity can 
be determined since all resonators receive the same flux^. 
After calibrating for the responsivity, we denote the "pri- 
mary channel" as the resonator showing the largest re- 
sponse, and the "opposite channels" as the group of res- 
onators showing the smallest response. 

As shown in Fig. 2, the pulse shape varies significantly 
between the primary and opposite channels. We find 
that the measured pulses, p(t), are well described by 
the sum of two components plus a noise term: p(t) = 
ApSp(t) + AoSoit) + n(t), where A^ is the amplitude of 
the prompt component that decays with the quasiparticle 
lifetime, Aq is the amplitude of the delayed component 
that decays with the phonon lifetime, and n(t) is a noise 



realization. The pulse shape for the prompt component 
and delayed component is given by Sp and 5o, respec- 
tively. The standard optimal filter formalism then gives 
the amplitude estimates for each pulse component for the 
ith channel as: 
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where the sum runs over the length of the trace, Sp^o and 
Pi denote the Fourier transforms of Sp^o and p^, J is the 
power spectral density of the noise and: 

M={ (3) 

The energy of the pulse is then reconstructed from the 
optimal filter amplitudes diS E = Eq X]i(^p,i +'^^o,i)/^i7 
where the latter component is weighted by the ratio of 
the template integrals, w = "^Sof^Sp, and the index i 
runs over the resonators. Since the two pulse components 
measure the quasiparticle creation rate convolved with 
the quasiparticle decay time, w gives the ratio of the to- 
tal number of quasiparticles created in these components, 
noting that the linear convolution does not affect the ra- 
tio. The relative responsivity of each resonator is given 
by Vi while the overall scaling, Eq^ is determined by cali- 
bration with x-ray lines of known energy. Figure 3 shows 
the reconstructed energy spectrum for phonon-mediated 
events from an ^^^I source. 

The frequency and dissipation response of the res- 
onators with temperature was used to determine a = 
0.075 ± 0.005 and A = 204 ± 21 /ieV^^'^^ Compar- 
ing the response versus temperature to the response 
from the source indicates that r]ph = 0.070 ± 0.011. 
Measurements of the frequency noise confirm that it 
is consistent with HEMT noise for T/v = 5 K. The 
quasiparticle lifetime, Tqp = 12.9 ± 1.2 /is, was deter- 
mined from the fall time of the prompt component of 
the primary channel pulses. Using these parameters 
and Eq. 1, the expected resolution for this device was 
(JE=iOA8 keV)^/r]read/Pt^ consistent with the measured 
baseline resolution for ^Jr]readlvt ^0.8. The phonon life- 
time was determined from the fall time of the opposite 
channel pulses to be ^50 /is, in agreement with the 1 ms 
fall times seen in TES-based detectors^ after accounting 
for the differences in substrate geometry and sound speed 
between the devices. 

Although the baseline resolution is consistent with ex- 
pectations, the measured energy resolution at 30 keV is 
cr^; ~ 1 keV, indicating that systematic variations in the 
pulse shape or amplitude with position are dominating 
the resolution. In addition, non-Gaussian tails of events 
lie between the expected spectral peaks. To remove 
poorly collected events and account for these systematic 
variations, a simple position estimator was constructed 
from the relative partitioning of energy between sensors 
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FIG. 3. (Color online) Observed spectrum from an ^^^I 
source. The reconstructed energy before (light histogram) 
and after (dark histogram) restricting to events interacting in 
the central portion of the substrate and applying the position- 
based correction to the energy are shown. The spectrum is 
fit to the observed lines at 29.5 (29.8%), 29.8 (53.1%), 33.6 
(10.2%), 34.4 (2.2%), and 39.6 keV (4.6%), where the num- 
bers in parentheses denote the expected absorbed intensities 
in 1 mm of Si. The thick line shows the fit to the total spec- 
trum, with best fit resolution of ge — 0.55 keV at 30 keV. 
The location of the 39.6 keV peak indicates there is a ^2% 
non-linearity over the 30-40 keV range, requiring a non-linear 
energy calibration similar to existing TES-based athermal 
phonon mediated detectors^^. (inset) Fit to the reconstructed 
energy spectrum for randomly triggered noise traces giving a 
baseline resolution of cr^; = 0.38 keV. 



for each event^'^^. For i resonators at locations {xi^yi) 
with energy Ei detected in each resonator, the "X energy 
partition" is defined as P^ = J2i ^i^i/ Z^i ^i^ and corre- 
spondingly for Y. Figure 4 shows the reconstructed event 
location using this position estimator, which has higher 
signal-to-noise than the corresponding estimator based 
on timing delays at 30 keV. Although further calibrations 
with collimated sources will be necessary to precisely de- 
termine the position resolution of the device, a rough 
estimate can be obtained from the overlap of the par- 
tition distribution for events with neighboring primary 
channels. The widths of these overlap regions indicate a 
typical position resolution of ~0.8 mm at 30 keV. 

For interactions occurring near the edge of the sub- 
strate, more phonon energy can be lost to the detector 
housing than for interactions in the center. Figure 4 
shows a position-based selection of events that removes 
edge events with poor collection, eliminating nearly all 
of the events in the non-Gaussian tails between peaks. 

In addition, a position-dependent correction for varia- 
tions in the reconstructed energy was applied, similar to 
the correction used in existing TES-based detectors^'^^. 
For each event with energy partition location (Px^Py), 
the set of n = 400 events with the closest partition lo- 
cation were selected. The spectrum for these "nearest 
neighbor" events was fit to determine the location of the 
29.8 keV peak, and the corrected energy for the event was 
determined as £;^^^^ = (29.8 keV/£;^9^8)^. where E^^^^ 
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FIG. 4. (Color online) Reconstruction of the interaction lo- 
cation from the energy partition. The coloring denotes the 
primary channel for each pulse. The black lines indicate the 
selection of events interacting in the center of the substrate. 
For comparison, the device geometry from Fig. 1 is overlaid. 



is the location of the peak maximum for the event's near- 
est neighbors. The resulting spectrum is shown in Fig. 3. 
The best-fit resolution is cr£;(30 keV) =0.55 keV. This 
resolution is nearly a factor of 2 better than the uncor- 
rected resolution for all events and is within 40% of the 
baseline resolution of aE{0 keV) = 0.38 keV. 

We are currently scaling this design to 0.25 kg sub- 
strates with Asub ^ 100 cm^. Existing TES-based 
phonon-mediated detectors of this size have typical en- 
ergy resolutions^ a^; = 1.0 keV at 20 keV, indicating that 
the devices presented above can already provide compa- 
rable energy resolution. The energy resolution of MKID- 
based devices can also be significantly improved. Increas- 
ing the phonon collection efficiency to r]ph = 30%, as 
obtained for TES-based designs with comparable metal 
coverage^, would improve the resolution by a factor of 
4. This collection efficiency was found to scale linearly 
with the film thickness and may be improved by reduc- 
ing phonon losses to the detector housing. Using res- 
onator materials with higher kinetic inductance or lower 
gap could improve the resolution, provided uniform res- 
onators can be fabricated. Tqp is also a factor of '^100 
smaller than results reported for 40 nm thick Al films at 
similar readout powers^^, indicating that a non-thermal 
quasiparticle population (possibly due, e.g., to stray 
light or particle interactions in the substrate) may be 
present in addition to readout generated quasiparticles. 
If so, better shielding of the device from external radia- 
tion could improve these lifetimes significantly. Finally, 
the development of lower-noise, broadband amplifiers'^, 
could provide up to an additional factor of 3 improve- 
ment in resolution. At the same time, MKIDs would 
provide less complex detector fabrication and higher res- 
olution position information to enhance background re- 
jection, simplifying the extension of these designs to the 
large target masses needed for future rare-event searches. 
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